%% 做出一个随机旋转角度的连接圆 

% function [输出变量] = 函数名称(输入变量）
function [x2,y2,locm,circle] = creatanconnectedCircle(circle_radius,x1,y1)

[xax, yax] = meshgrid(1:227, 1:227);
rax = xax+1i*yax;

loc1 = x1+1i*y1;

% 圆2位置
rt1 = abs(rax-loc1);%rt1<circle_radius就是圆1;下行命令是rt1<4*r，就是为了找出圆2运动轨迹围成的大圆，然后再用"bwboundaries"抠出轨迹
loc2ind = rt1<4*circle_radius; % 第2个圆的圆心运动轨迹的大圆
loc2bound = bwboundaries(loc2ind); % 第2个圆的圆心运动的轨迹（loc2ind的圆周长），是cell
loc2bound = loc2bound{1};% 将cell提取出来，形成矩阵，方便下面提取
loc2 = loc2bound(randi(size(loc2bound,1),1),1:2); % 提取圆2的圆心运动坐标的范围，随机挑一个坐标作为圆2的坐标
loc2 = loc2(2) + 1i*loc2(1);% 实部对应x（列），虚部对应y（行）

% 连接圆的坐标参数
x2 = real(loc2); % 圆心2的横坐标
y2 = imag(loc2); % 圆心2的纵坐标

%中点位置
xm = (x1+x2)/2;
ym = (y1+y2)/2;
locm = xm + 1i*ym;

circle = [x1,x2;y1,y2];

end